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I. INTRODUCTION 



The nature of QCD string formation between static sources was studied on the 
lattice I^J^ and analytically Q. From these investigations it was shown that the 
string consists of a predominantly color-electric longitudinal field. At the critical 
temperature Tc this electric field disappears and and at the same time the deconfined 
phase with color-magnetic condensate sets in. This effect was predicted theoret- 
ically in 1^ and also seen in lattice measurements [Q. At the same temperature 
Chiral Symmetry Breaking (CSB) for light quarks is found to disappear jsj, which 
indicates that there is an intimate connection between the string formation and 
CSB. In the case of heavy quark systems CSB occurs due to the quark mass and 
confinement can be described as the area law of the Wilson loop. How CSB and 
confinement are explicitly realized for the light quark system and what equation 
describes its dynamics is an interesting and open problem. 

It is the purpose of the present paper to study this issue in the simplest dynamical 
example - in the system of one light quark and a heavy antiquark. This allows us to 
describe the dynamics of light quark (its propagator) in a gauge-invariant manner, 
while physically the light quark is expected to be confined at the end of a string 
connected to the static source. Applying the formalism of field correlators (FC) 
I^J^, we derive the effective quark Lagrangian, containing any number of quark 
operators multiplied by field correlators. To proceed further one can use the limit 
of large Nc and write down the Dyson- Schwinger equations with the mass operator 
expressed through the Green's function. The resulting equations arc nonlocal and 
nonlinear. It is not clear from the beginning how confinement and CSB would 
manifest themselves in the solution of these equations. Some hint was provided 
m Ref. H using a relativistic WKB analysis |9j, where it was shown that at large 
distances from the heavy source the dynamics of the light quark is described by the 
Dirac equation with a scalar linear confining interaction, which leads to CSB. 

In this paper we examine the properties of the Green's function of the color 
singlet qQ system, where the antiquark is treated in the static limit. In section II 
we formulate the full form of nonlocal and nonlinear equations for the light quark 
propagator and its eigenfunctions and study the behaviour at all distances. We 
also take into account both perturbative and nonperturbative contributions to the 
interaction kernel. As a result our equations contain both a confining interaction 
and color Coulomb part. Similar to the heavy quark situation we argue that the 
string formation for low angular momentum is of a color-electric nature. Moreover, 
the confinement of the light quark to the heavy one is shown to be of Lorentz scalar 
type. In section III the resulting nonlinear equations are studied numerically. The 
energy spectrum and the structure of the low lying eigenfunctions are presented. 
We in particular study the B and D meson spectrum. We compare them with 
experimental data for B, D mesons and results of other calculations, exploiting for 
this purpose the expansion of Heavy Quark Effective Theory (HQET). In section 
IV the chiral condensate of the light quarks {qq) is determined by taking the limit 
of the Green's function S'(x, y), with both x and y tending to zero. In this limit the 
heavy quark is turned off and the condensate {qq) ~ 5(0, 0) should have a value not 
depending on the presence of heavy quark. 



II. DERIVATION OF EQUATIONS 
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A. Dyson- Schwinger equations 



In this section we give an outline of the procedure to obtain the Dyson-Schwinger 
equations for the color white qQ system. Our starting point is the gauge-invariant 
light quark Green's function S{x, y) in the presence of a static heavy antiquark 
placed in the origin. In the static limit, the heavy quark can be treated as an 
external source. Assuming the Euclidean metric and letting T = {x — y)4, the heavy 
antiquark propagator in the modified Fock-Schwinger gauge pO[ | is proportional to 
the parallel transporter, namely 



Sq{A) = h{x, y) Pexpiig [ A^ir = 0, T)dT) 

Jo 



(1) 



with 



/i(x,y) = -^(3)(x- y)[(l + 74) e-^^^eiT) + (T ^ -^,74 ^ -74)]. 

This acts as a static source situated at the origin. As a result the proper limit for 
toq — > 00 of the Green's function of the qQ system can be defined as 

S{x,y) = mx)Uix,y)^Jj+iy)). (2) 

Here we have to average over the gluon fields A and light quark fields ip, while 11 
contains the parallel transporters between the end points 

n(x, y) = 0(x, X4; 0, X4) 0(0, X4; 0, 2/4) 0(0, 2/4; y, Va) 

with 



(x, a;4 ; y , y4) = f exp ig ( A_ 

J X 



The averaging over the gluon fields A has to be done over the perturbative and 
nonperturbative gluon fields and contributions respectively, where the to- 
tal gluonic field = i?^ + a^. We now apply the method of field correlators 
(FC), which was developed in a series of papers [^,0 to derive the effective quark 
Lagrangian from QCD. Let us first consider the effects of averaging over nonper- 
turbative (NP) field Bf^. We may write for the partition function 

(3) 

with j{n) = jp,„{xn) — i'^ {xn)lfinip{^n) and B{n) — _B^^(x„). To write the corre- 
lator {{B...B)) for the gauge-invariant situation corresponding to the color white qQ 
system one can use the modified Fock-Schwinger gauge to express the correlator 
of Bix) through FC: 

7V(1, ...n) = {{B{l)...B{n))) ^ [ dx{l)... dx{n){{F {!)... F{n))). (4) 



The effective interaction kernel in Eq. (|^) can now be used to write a Dyson- 
Schwinger-type equation for the quark Green's function S. To simplify matter, 
one can consider the large Nc limit, in which case the connected self-energy kernel 
M{x,y) is obtained from Eq. (|^) by replacing any pair of adjacent ?/;-operators by 
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'^aa{x)^^p{y) NcSai3{x, y), (5) 

where aa are color and Lorentz index respectively (for details of derivation see 1^). 
As a result one obtains the equation for S{x,y) 

{-icjlx - im)S{x,y) ~ i / M{x, z)S{z,y)d'^z ^ d^^\x - y), (6) 



where the kernel M is expressed through N as 

tM{x, y)= N^l\x, yh^S{x, yh, + (7) 

OC „ 

^ / d^^X2..J'^x„^lJf,^S{x,X2)-f^,^...J^,„^^S{x„-l,y)J^,,^N^^'^l^^^{x,X2■■■Xn-l,y). 

n=3 

The system of equations (6-7) is exact in the large Nc limit and is well defined 
provided all NP correlators {{F {!)... F{n))) are known. 

Evidence has been found in recent accurate lattice calculations [0 of static po- 
tentials in different SU{3) representations, that the contributions of the higher 
correlators {{F {!)... F{n))) for n > 2 to the planar Wilson loop are small. In partic- 
ular, these terms are found to contribute only around a few percent of the dominant 
Gaussian correlator |l2|. Hence, the Gaussian Stochastic Model, based on the low- 
est correlator ((F(1)F(2))) is expected to be a good approximation. In view of this 
we assume in this paper that the Gaussian approximation holds and we keep simi- 
lar as in Ref. |^ only the first term in Eq. (^). One can parametrize the Gaussian 
correlator according to Q as 

((F(l)^^(2)» - ■^Jr{F^x{x)^{x,0)F,40)^{0,x)) ^ Dix){S^.Sxa - S^^S.x) + 

(8) 

where only D{x) is responsible for confinement and it contributes to string tension 
a, while A^^^ is a full derivative. As a consequence, the latter contributes to the 
perimeter of Wilson loop and (f>{x, 0) ~ P exp ig B^dz^. From this we find, that 

the n'i^J can explicitly be written in the gauge (lO| as 

N^j^J{x,y) ^ {Sf_,^Sik - d^kS,yi) duia^^{u) j dvka^{v)D{u - v), (9) 







where 04 (it) = l,ai{u) — 

In Eq. (|p) only the nonperturbative confining piece of the Gaussian correlator (H) 
is retained, since the perturbative part and A*-^) do not produce neither the string 
(confinement) nor CSB In the case of a nonrotating string the terms in Eq. (7) 
with space components 7^ are suppressed by powers of velocity of the endpoints 
of the string. In what follows we shall keep for simplicity only the component 
N^^' = J{x,y) of Nf^J . Hence the kernel M is proportional to the FG of the 
color-electric field Ei = Fa. It is the dominant part of the string. Golor-magnetic 
components are neglected in this first step and can be considered as a correction. 

In this way one arrives at the system of equations |^ where we keep the same 
notation M{x,y) for the retained piece of the kernel 

iM{x,y) ^ J{x,y)'^4S{x,y)j4 (10) 
i-i^x -tm)Six,y)^i J M{x,z)S{z,y)d^z = S'^^\x-y). (11) 
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B. Partial wave reduced equations 

The kernel J{x,y) in Eq. ( |l0| ) depends on the time as ^"^^^ . Therefore in the 
Umit of Tg — !■ it becomes local in time. As a result, the Fourier-transform of M 
in the fourth coordinate does not depend on the momentum p4 for Tg — s- 0. The 
corresponding eigenfunctions ^„(x,p4) = V'n(x) and eigenvalues £„ = SniPi) of the 
l ight q uark in the white light-heavy configuration can readily be obtained from Eqs. 



(10 111) in the discussed instantaneous limit. We find 
c 

t dx 



{--^ + (3m)^IJnix) + P I M(x,z)V„(z)d^z = £„7^„(x) (12) 



Af(x, z) ^ J(x, z)/3A(x,z), A(x,z) = ipk{x)signek'4)'l {z) (13) 

k 

The spherical-spinor decomposition of "i/'n for the total and orbital angular momen- 
tum channel j, / = J ± ^, 

yields equations for the partial waves 

^ - -F„ + (£„ - m)G„ - MnG„ - iKh^Fn = (15) 
ar r 

dG 

+ -Gn - {en + m)Fn - M22Fn + lAhlGn = (16) 

ar r 

with K = ±(j + i). Clearly the kernels Mik are nonlocal in space, i.e. 

/•OO 

M,fcG„=/ Mll{ry)Gn{r')rr'dr' (17) 



with Mil = {^.M\Mik\^:ji^M) and M,k is given by Eq. (|T|). 
Equations dTa-nW) are invariant under the transformation 



-e„, K ^ -K, G„ En (18) 



which also yields Mn <^ M22,Mi2 ^ —A'hi- The symmetry ( |1§| ) implies that the 
spectrum is symmetric in £„ <-> — e„, which is a property of a scalar interaction. 
The Lorentz scalar nature of the confining interaction has the nice feature that it 
doesnot lead to instability problems in the Dirac equation [Q. Moreover, nontrivial 



solutions of Eqs. (|15| 
We solve Eqs. (^ 



-16), if they exist, signify spontaneous CSB. 
^16), using the relativistic WKB approximation for the kernel 
Mifc. To simplify the calculations the Gaussian form for D{x) was used (since all 
observables are integrals of FC, its explicit form is not essential at large distances, 
provided the FC have a finite range Tg and it yields the same value of the string 
tension cr) 

D{u) = D{0)exp{-uy4T^), D{0) ^ ^ (19) 

with a = 0.2 GeV^, and Tg — 0.25 fm, taken in accordance with the lattice 
measurements |l|,[l3[. 
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As the reference basis we take the WKB sohitions of the Dirac equation for 
the local linear confining potential ar, and the WKB computed kernels M and 



A(x,y) in Eq. (13) are determined by explicit summation over eigenstates. It was 
checked by an independent calculation ||l^ that the relativistic WKB procedure 
yields eigenvalues of the linear potential with accuracy better than one percent. 

The general structure of M and M can be derived from Eqs. (|l^-|lj). One can 
write M as a 4 X 4 matrix as follows jl^ 

M = M^°h + M^'^a, + M('*)74 + M«7„ (20) 

where 7^, 74 are usual Dirac matrices, i = 1,2,3 and Ui — ^ q * |^ 

The same representation holds for the WKB approximated M . From the WKB 
analysis we find that Af'"' is the only growing kernel. It behaves asymptoti- 

cally as 

M(")(x,y)=a^i±^5(3)(^_y) (21) 

where ^'^•^^(r) is a smeared 5 - fmiction with the range of nonlocality decreasing 
asymptotically with growing |x|, |y|. The term M^'^ is proportional to the angular 
momentum L and asymptotically it behaves as 0{\/x). One can also prove that 
M'^^\m\"^ do not grow at large x,y |l^. Hence in all problems where large dis- 
tances are dominant one can consider only the first term in Eq. ([20|). Using this 
approximation we get for the kernel ( p^ 

M = M(") (x, y)I = J(x, y) A(x, y)I (22) 

with 

J(x,y) =a^/(x,y), /(x,y)== / ds dte (23) 

Jo Jo 

and where A is now a scalar quantity. 

A curious feature of the considered equations is the way how the string connecting 
the light quark to the source is being created. Actually, if one takes only lowest 
partial waves inside the kernel M (i.e. in A(x, y), Eq. dl^)), then the effective 
potential in Eqs. ( p"5[]l6| ) is not confining. If one however sums up over all angular 
states and radial excitations in A, then the resulting A(x, y) is a smeared S- function 
leading to the quasilocal confining kernel M. E.g. using eigenfunctions for the local 
case, A can be computed quasiclassically to be 



A(x, y) = — ^ ^-====== , xy = xy cos 9 (24) 

Ztt-^ ^[x - yY + e^xy 

In Eq. ( p^ ) one can clearly see that K{x,y) is a normalized smeared (5-function, 
with smearing radius in |a: — j/| being ^J^ . For large distances it is nonvanishing 



only in the forward direction. 

Insertion of this A into M, Eq. (^2|), produces linear confinement due to the kernel 
A(x, y), as is given by Eq. (P4| ) (while simply averaging the contribution from each 
individual orbital in Eq. (|l3| ) over the angle between x and y would produce no 
confinement at all). The computed kernel M{x,y) turns out to be nonlocal, but 
very close to the linear potential at large distances. Indeed, the effective localized 
potential defined as 
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Veffir) = / M^"\r,x)rxdx 



(25) 



approaches at very large distances, i.e. for cr^^^r > 200, a linear dependence with 
a slope given by a. However at shorter distances Ve// looks also linear over a 
relatively large region with a slope almost the same as a, reflecting the presence of 
a small local curvature. In particular, we find that in the region 5 < cr^/^r < 20 the 
effective potential can reasonable well be described by Va{r) — O.Qar — 1.80-1/2. In 
Fig. 1 are shown the results up to r = 20 in units of cr^/^. 

Since it is only the higher states and large distances which are important in the 
creation of this S - function-type behaviour of A(a;, y), and since the WKB method 
does well for high states and at large distances, one can clearly conclude, that linear 
confinement should be obtained if one sums over all exact solutions of Eqs. ( p"5[ - |l^ ). 
Hence this should be a property of the exact solution. 

The property, that the kernel has a focussing effect in the forward direction can 
be used to get a somewhat simpler form. For this purpose we may also use Q 

A(x, y) = Ko{a)S{cose^ - cosdy), a = (J^/xy\x - y\ (26) 

TT^'^/xy 

The eigenvalues and eigenfunctions in Eqs . (|l 5^161 ) have been determined using 
the kernel M(x,y), given by Eqs. ( p^ ) and (|26|). Some results are shown in Table 
1 and Fig. 2. In our present study we have taken m = 0. Note, that M{x,y) is 
approaching a local linear potential at large distances, x,y > cr^^/^, which justifies a 
posteriori our choice of the reference basis. Due to the nonlocality of the interaction 
the predicted spectrum is found to be different from that of the linear potential, 
valid at large distances. Moreover, comparing the level structures of the J = 
channel as obtained using the kernels (^4|) and (^) we see from Table I, that they 
are qualitatively very similar, corroborating that there is indeed a strong forward 
focussing effect in the quark propagator. From Fig. 2 we see that for all L values 
the higher radially excited levels are close to the predictions of linear potential 
Vo{r) = 0.9ar — 1.8a^^^, in agreement with the fact that the interaction at large 
distance can indeed be described by a local linear potential. On the other hand the 
nonlocal kernel predictions for the low lying states clearly deviates strongly from 
those of the (shifted) linear potential. Hence the nonlocal nature of the force does 
affect the spectrum in an essential way. 

The eigenfunctions for the nonlocal kernels look qualitatitively similar to the 
corresponding ones of the shifted linear potential. In Fig. 3 are shown the ground 
state and first excited state for the J — channels. Altough the differences are 
substantial for these low lying states, the agreement for higher excited states is 
considerably better. Moreover, we find that the large distances and high states of 
the WKB states agree well with the corresponding eigenfunctions. 



C. Inclusion of perturbative exchanges 

Till now we have considered only NP part of the gluonic field, S^. In this section 
we include the perturbative part, a^, and neglect for simplicity the interference 
terms. Therefore the effect of is accounted for in the appearance of an additional 
factor in the partition function (|^), namely 

Zj — Z^NP'^pert, Zp^rt — •> •> )a, I ) 

where the second term in the exponent of ( p7[ ) corresponds to the interaction of 
the perturbative part of the gluon field with the static antiquark. We have used in 
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(27), that due to the 'tHooft identity ||17[] one can average independently over i?^ 
and a^. 

The result of averaging yields a new additive term in Lg/ / , Eq. , 

Lc = .g j dxi^+{x)A^''\x)il}{x), (28) 

where we have defined 

4-)(a;) = -ig j dz^ < a4(x)a4(^4) >= 6^4 ^'^^^'^ . (29) 

The presence of Lc in Eq. (|^) does not influence the derivation of basic Eqs. ( p^fjll] ). 
The only difference is that Eq. (O) assumes the form 



{-id ~ gA'ix) - im)S{x, y)~i M{x, z)S{z, y)d^z = S^^^ (x-y). (30) 



Eq. ( p^ does not change and the kernel J{x, y) contains as before only nonper- 
turbative contributions. Note however that S{x,y) in M{x,y) in Eq. ( |l0| ) now 
contains also perturbative gluon exchanges. This is a new type of interference of 
perturbative and NP terms, which appears irrespectively of our neglect of this in- 
terference within the averaging procedure over and a^. In other words another 
class of diagrams is responsible for this interference. 

Correspondingly in the static equations (na) one should replace 



(3m j3ra ^—^ (31) 

|z| 

The equations for the partial waves ( [l5| - L6| ) are modified due to the presence of the 



color Coulomb potential V{r) in a simp. 



e way. Since 



y(,)=_^ (32) 

is local and a Lorentz vector, it always appears in the combination e„ — y(r). Hence 
one has instead of Eqs. (|l5|-[16|) 

^ - -F, + (e, - V{r) - m)G, - M^G - iMi^F = 0, (33) 
dr r 

-T^ + -G,- (£. - V{r) + m)F, ~ M22F + lAhiG = 0, (34) 
dr r 



where we have denoted 
G 



F 



<i^mkW'> (^fj^^^^ywdw. (35) 

Here Mik is defined as in Eq. (^3|) and the matrix A^^ in Eq. ( p^ involves the sum 
over all states, including positive and negative £„. There in section IIB we have 
exploited the symmetry (|l^). However Eqs. (|3^p4[) are invariant under another 
transformation, namely 

£n ^ -£n, V{r) ^ - V{r), n ^ n, G„ Fn- (36) 

Now the sum over negative £„ can be expressed through the corresponding sum 
over positive £„ with exchange Gn <-> Fn as before, but also with the inversion of 
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sign of Coulomb interaction, i.e. Coulomb attraction for positive £„ is replaced by 
Coulomb repulsion for negative e„. 

In what follows we shall denote wave functions of the positive energy states with 
repulsive Coulomb with the sign of tilde: G^^Fy. Then using (13) the matrix /JA^fc 
can be written as a sum over only positive follows 

^ . pm' ^ J_ ( G^,G*^, - F^,F*,, -i{G^,F*, - Ff,G*^,)\ 

' ^^.i;r>oV-*(^''^M'"G,i?;')' g,g;,-f,f;, )■ ^ > 

Since fHA is exactly the combination which enters the mass matrix ( [Tsl ) , one can list 
in ( |37| ) scalar and vector (proportional to (3) parts: 

M = Msl + M^P + AM, (38) 

where AM contains spin-dependent terms, which can be considered as in section 
IIB, while Ms, My are 

M,,,.^G J2 [G,G;-F^F;±iG,G^-F^F;)] (39) 

where 

C=iv^T,i^(0)^/(x,y), f{^,y)^ j\sj\texpiy-^^^^^^y (40) 

From (^) it is clear that the vector part M„ is only due to the presence of Coulomb 
interaction. Corrections at large distances due to the vector part can be treated 
again in the relativistic WKB. A rough estimate of My at large r yields 



My _ a 
Ms ^ ar 



(41) 



and hence can be neglected at large enough r. 



III. NUMERICAL SOLUTIONS OF EQUATIONS AND COMPARISON 

TO B, D MESONS 

We have performed numerical studies of Eqs. ( ^SfjS^ ) with the kernel (^^ for 
different values of the quark mass m and different values of Tg. To simplify 
calculations only the dominant part of the mass operator Mik was retained, i.e. 
Mil = M22 = M^'^\ while Afi2, A/21 have been neglected. For Af'"' the representa- 
tion ( p^ ) was used 

A/(o)(x,y) = J(x,y)A(x,y)/, 

where A is taken to be the kernel (|2j). Results of our calculations for the ground 
state energy are listed in Table 2. One can see a rather sharp change of energy 
when Us changes from to 0.3 and when Tg is changing from to 0.25, while 
further increase of as or Tg does not produce such a strong dependence. 

Solutions of our equations ( ^3[]3^ ) can be compared with physical states oi B, D 
and Bs, Ds mesons. To this end one should have in mind that in Eqs. (|33| - ^4|) the 
static approximation for the heavy quark b, c was used, and hence all corrections 
0(1/771^) with n> 1 are neglected. 
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One can exploit at this point the HQET expansion for the mass mn of heavy-hght 
boson |l|,|l9| 

= ™q(1 + — + TT-ri^i + + 0(l/m3 ), (42) 

where A„ are free parameters, depending on dynamics, and dn is the hyperfine 
spHtting parameter. It is clear from the preceding that eigenvalues of Eqs. (33-3^) 
yield the value A, which depends on the quantum numbers of the state, 

A(j,;,n^) =£„(j,0 (43) 

Consider now the results of the present approach, i.e. solutions of Dirac-type 
equations ( ^3[]3^ ). In the local case (Tg 0) when the kernel M reduces to the 
linear potential trr, we have 



and 



a(^°=) = 0.690 GeV («« =0, ct = 0.18 GeV^) (44) 



A^'''^^ 0.493 GeV^ (a^ = 0.3, a = 0.18 GeV^). (45) 



This should be compared to the nonlocal case 



and 



A^°"'°"^ = 0.288 GeF (a, = 0.3, o- = 0.18 GeF^). (47) 



These latter values arc in general agreement with the results of the QCD heavy- 
flavour sum rules p0|j2l[ | 

A = 0.57±0.07 GeF (48) 
and more recent analysis from semileptonic B decays |2^ 

A = 0.39±0.11 GeV^ (49) 

Another interesting comparison is with the experimental values of the i?-meson 
mass (the term A2 in Eq. ( [42| ) can be determined from the B* — B mass difference). 
Using Eq. (^) and Ms = ^Mb,+Mb ^ 5 3;^2 GeV one can estimate (neglecting 
Ai) the pole mass of the 5-quark to be mi,{pole) = 5.0 GeV, which is in reasonable 
agreement with the analysis of the quarkonium spectra in |23[ . 

A similar analysis can be done for the Bs meson; the corresponding values for A^ 
with TTis = 0.15 and 0.20 GeV, are A^ - A = 0.084 and A^ - A = 0.115 for = 0.3 . 
One can compare these values with the mass difference Bs,B'-\ AMs{B) = (0.090 ± 
0.0038) GeV. These numbers for A can be compared with those in Table 3, where 
also results of lattice calculations and of the constituent quark model (CQM) 
are given. 
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IV. CHIRAL CONDENSATE 



As a check of CSB in our Eqs. ( |l5| - |lq ) we have computed the chiral condensate, 
which can be expressed through the eigenfunctions as in |^ (to simphfy matter we 
disregard in this section perturbative contributions). 

N °° 

i^i) = -^j2^iA-r-(Bn)% (50) 

n=0 

where — { '^"^^ )r=G, = (■^^^7^)r=o, and Gn,Fn refer to solutions with k = 
— 1,1 = and K = +1,1 = 1 respectively. In the local linear potential case the 
values of A~ , have been computed in the WKB method |^ and shown to yield 
a monotonically divergent series {qq) — J2n ■ 

It can be argued (using Eqs. (8-9) from |^), that the nonlocality of the kernel 
M in space-time, present by definition in Eq. (8) improves the convergence of the 
series and yields a finite result for (qg) . We have found A'^ , 5+ from the solutions 
of the nonlocal equations ([l5|-p^ with the kernel M and compared them with the 
local case, when M reduces to the local linear potential. Results are shown in Table 
4. 

One can see from the results, that in the nonlocal case the magnitude of s„ = 
(A~)^ — (i?+)^, is clearly diminished as compared to the reference local case, and 
is of reasonable order of magnitude. From the obtained sequence of s„ we get that 
(qq) = —0.5(7^/^ and — 0.7fT^/^ in the nonlocal cases of the kernels (21) and (22) 
respectively. Adopting a value of cr = 0.2 GeV"^ we find {qq) = —(350 MeVY 
and —(400 MeV)^ respectively, to be compared with the usually acceptable value 
of —(250 MeV)^. However convergence is still slow as seen from Table 4 and 
the converged values are somewhat higher. We have checked that convergence is 
somewhat improved when one takes into account the intrinsic nonlocality of the 
kernel M in x, y. To this end we have modified the kernel M obtained from WKB 
analysis, replacing S in Eq. (21) by a Gaussian factor 

Nexp{- ^^~J^^ )6(cos0r - cos0y) (51) 

and studied the sequences of s„ as functions of the nonlocality range a. Results 
are shown in Table 4 for two values of a = 0.3cr^/^ and O.ba^^^. The strength 
N is chosen such that numerically the slope of ct = 0.2 GeV^ is reproduced for 
large distances. The condensate values varies in the considered region substantially, 
showing that effects of lonlocality are important. 

The slope of the effective potential Veff{r), determined by N, strongly depends 
on a for a = Tg. We believe that the reason for this lies in the fact, that the 
chiral condensate (qq) depends crucially on the nonlocality both in time components 
of Af (x, y; X4, j/4) and in spacial components. The first nonlocality was however 
disregarded in Eqs. (|l^|l3|) , when the p4 - dependence was omitted in M and ipn 
(the static limit). It was indeed shown in Ref. that taking this dependence into 
account significantly improves convergence of the sum in Eq. (|5^) . The full account 
of this effect requires solution of time-dependent Eqs. ( p^pT| ), which is numerically 
a much more difhcult problem. 



V. CONCLUSION 

We have studied the confining and CSB properties in the system of one light 
quark and one static antiquark. The effective mass operator is written explicitly for 
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large Nc, as a sum over vacuum field correlators. Keeping only the Gaussian field 
correlator, we have obtained a closed system of equations in the limit of large Nc- 
Our results support the presence of a Lorentz scalar linear confinement for the light 
quark, which signifies CSB for this system, and yield eigenfunctions and eigenvalues 
for the heavy-light system containing both confinement and CSB. 

As a direct evidence of CSB we have computed the chiral condensate, which 
appears to be of the correct sign and having the proper large Nc dependence. Our 
result yields a reasonable order of magnitude of (qq), provided convergence of the 
sum is achieved. At this point it is useful to compare the CSB picture of the NJL 
model and our approach. In the NJL model confinement and string are absent 
and CSB may occur due to the condensation of qq pairs in the scalar channel. In 
our case, being the large Nc approximation of the real QCD, a string is built up 
between light and heavy quark, which depends not only on light quark coordinates 
X, y, but also on the distance from them to the heavy antiquark. In the presence 
of confinement, the phenomenon of CSB is due to the spontaneous creation of the 
scalar string, which is forbidden by chiral symmetry. 

Eigenvalues e„ and eigenfunctions obtained numerically for lowest states, repre- 
sent the leading contributions of the HQET expansion in powers of l/mg. Results 
for the energies e„ in our method are compared of the lattice and QCD sum rule 
calculations, and also with experimental extraction of e„ = A(n), showing an overall 
agreement with the B and D meson masses. 
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FIG. 2. Level Structure calculated with the Dirac equation using the kernel ( pq ) (solid 
line) in comparison with the predictions of the linear potential Vo{r) = 0.9r — 1.8 (dashed 
line). States carry the quantum numbers Lj, with L, J being the orbital and total angular 
momentum. 
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FIG. 3. Eigenfunctions of the ground state and the first excited state for the Si and 
Pi channels. The solutions correspond to using the kernel ( p6| ) (solid line) and the shifted 
linear potential Vo(7') — 0.9r — 1.8 (dashed line). 
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Table 1. Energy eigenvalues for Eqs. (13-14) with J = ^ and the kernels, given 
by Eq. ( p^ ) and Eq. (p6|), as compared with those of the Dirac equation for the 
linear local potential V/i„(r) = r. 



V ■ 



K = -l 



1.619 
2.603 
3.291 
3.855 
4.345 
4.784 
5.186 



K = 1 



2.294 
3.031 
3.626 
4.138 
4.594 
5.008 
5.334 



kernel (24) 



K = -1 



0.925 
1.719 
2.277 
2.740 
3.144 
3.507 
3.838 



K=l 



1.472 
2.056 
2.541 
2.964 
3.342 
3.685 
4.002 



kernel 



K = -l 



0.969 
1.765 
2.334 
2.809 
3.226 
3.603 
3.950 



K=l 



1.516 
2.113 
2.608 
3.042 
3.432 
3.790 
4.123 



Table 2. Ground state energy eigenvalue (in units of y/a) for Eqs. (33-p4|) with 
as — 0, 0.3 and 0.39 and quark masses m = 5 MeV, 0.15 GeV and 0.2 GeV (upper, 
middle and lower entry) for different values of Tg, Tg = 0, 0.25, 0.5 and 1 (in units 
of l/^/^). 



as 





0.25 


0.5 


1 





1.628 
1.886 
1.978 


0.985 
1.225 
1.314 


0.979 
1.217 
1.305 


0.907 
1.145 
1.233 


0.3 


1.163 
1.378 
1.456 


0.684 
0.884 
0.959 


0.679 
0.877 
0.951 


0.628 
0.826 
0.900 


0.39 


1.004 
1.201 
1.272 


0.585 
0.768 
0.837 


0.580 
0.761 
0.830 


0.536 
0.717 
0.786 
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Table 3. Energy eigenvalues A of the heavy-light system in the static heavy quark 
approximation obtained in different approaches. 



Refs. 


Method 


A (GeV) 


20 


QCD sum rules 


> 0.5 


21 


QCD sum rules 


0.4^0.5 


24 


Lattice 


0.18 ±0.03 


22 


expcrim. 


0.33 ±0.11 


25 


QCM 


0.35 


26 


QCM 


0.5 0.6 


27 


Rel. QCM 


0.386 


this work 


Nonlin. Dirac Eq. 


0.287 



Table 4. The difference s„ = — for n = 0, 1, ...6, in case of the nonlocal 
kernels (|2^) and (|6|) and corrected for a normalized Gaussian nonlocality (|l]) with 
a range of a =0.3 and a=0.5. For comparison the results are shown for a local linear 
potential Vun — r . 



n 




V = r 


kernel (24) 


kernel (26|) 


a=0.3 


a=0.5 





0.79 


0.42 


0.50 


0.15 


0.23 


1 


0.51 


0.21 


0.34 


0.04 


0.12 


2 


0.41 


0.12 


0.19 


0.02 


0.10 


3 


0.35 


0.10 


0.16 


0.02 


0.09 


4 


0.31 


0.08 


0.11 


0.01 


0.09 


5 


0.27 


0.07 


0.11 


0.01 


0.08 


6 


0.26 


0.06 


0.09 





0.07 
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